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Abstract.  From  the  analysis  of  VLBI  observational  data  compiled  by 
USNO  (U.S.  Naval  Observatory)  from  MJD  44089.994  to  51618.250  (Mc¬ 
Carthy,  2000),  we  showed  that  a  strong  peak  around  —400  sidereal  days 
in  the  spectrum  of  its  differences  from  the  IERS96  nutation  theory  could 
be  explained  by  adding  a  model  Free  Core  Nutation  (FCN)  term  in  the 
form  of  a  single  damped  oscillation.  Then  we  developed  a  new  analytical 
theory  of  the  nonrigid-Earth  nutation  including  the  derived  FCN  model. 
We  adopted  RDAN98  (Roosbeek  and  Dehant,  1998)  as  the  rigid  Earth  nu¬ 
tation  theory.  It  was  convolved  with  a  transfer  function  using  numerical 
convolution  in  the  time  domain  (Shirai  and  Fukushima,  2000).  The  form 
of  the  transfer  function  was  the  same  as  that  of  Herring  (1995).  How¬ 
ever,  its  free  parameters  such  as  the  complex  amplitude  and  frequency  of 
the  FCN  were  readjusted  by  fitting  to  the  above  VLBI  data.  Even  after 
truncating  the  forced  nutation  series  so  as  to  contain  only  180  terms,  the 
WRMS  (Weighted  Root  Mean  Square)  of  the  complex  residuals  for  the 
new  nutation  series  is  0.312  mas,  which  is  significantly  smaller  than  0.325 
mas,  that  of  the  IERS96  nutation  theory.  As  for  the  FCN  term,  we  esti¬ 
mated  its  oscillatory  period  as  —430.8  ±0.6  sidereal  days,  and  its  Q-value 
as  16200  ±  1600.  Also  we  estimated  the  correction  of  the  precession  con¬ 
stants  as  — 0.29297±0. 00047"/ cy  in  longitude  and  — 0.02430±0. 00019"/ cy 
in  obliquity,  respectively. 


1.  Introduction 

The  nutation  has  two  modes,  the  forced  and  the  free  modes.  The  forced  nuta¬ 
tion  is  due  to  the  lunisolar  and  planetary  tidal  forces.  The  free  mode  appears 
because  the  Earth  has  a  rotating  and  ellipsoidal  fluid  core.  Thus  this  mode  was 
named  as  the  free  core  nutation  (FCN)  by  Toomre  (1974).  However  the  current 
IAU  nutation  theory,  which  is  called  the  IAU1980  nutation  theory  (Seidelmann, 
1982),  does  not  include  the  FCN  term.  This  is  because,  at  the  time,  the  FCN 
had  not  yet  been  observed  directly  optically  (Sasao  and  Wahr,  1981). 
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The  IERS96  nutation  theory,  which  is  in  the  best  agreement  with  the  VLBI 
observation,  also  does  not  contain  the  FCN  term  (McCarthy,  1996).  See  Figure  1 
for  the  difference  between  the  VLBI  observation  and  IERS96  nutation  theory  for 
the  period  MJD  44089.994  to  51618.250.  There  remains  a  strong  peak  around 
—400  sidereal  days  in  its  spectrum,  shown  in  Figure  2.  This  is  probably  because 
the  strong  peak  in  the  spectrum  of  the  residuals  is  believed  to  be  explained  by 
considering  external  excitations  which  vary  with  the  periods  close  to  that  of  the 
FCN.  As  shown  by  Gegout  et  al.  (1998),  the  atmosphere  is  a  good  candidate  for 
such  an  exciting  force,  which  means  that  the  amplitude  of  such  peaks  cannot  be 
predicted.  However,  the  strong  peak  may  be  explained  appropriately  by  adding 
a  model  FCN  term.  In  this  paper,  we  will  try  to  improve  the  analytical  theory 
of  the  nonrigid-Earth  nutation  by  adding  a  simple  model  term  of  the  FCN. 


2.  Model  of  Free  Core  Nutation 

We  first  fixed  the  forced  nutation  theory  as  the  IERS96  nutation  theory.  We  set 
the  frequency  of  the  FCN  as  that  of  Herring  (1995).  As  for  the  model  of  FCN,  we 
considered  the  following  two:  (A)  no  excitation  during  the  observation  (model 
A),  and  (B)  frequent  excitations  during  the  observation  (model  B).  Model  A 
is  mathmatically  expressed  as  a  single  damped  oscillation,  where  the  damping 
is  due  to  the  viscosity  of  the  outer  fulid  core  and  the  magnetic  coupling  of 
mantle-core  boundary.  Model  B  is  regarded  as  a  simple  harmonic  oscillation 
without  damping  whose  (real)  frequency  u  is  close  to  the  real  part  of  the  FCN 
frequency.  In  both  models,  we  estimated  one  complex  amplitude  of  the  FCN  by 
weighted  least-squares  fitting  to  the  residuals.  In  model  B,  we  also  derived  the 
real  frequency  u  by  fitting.  In  both  models,  the  strong  peak  in  the  spectrum 
around  —400  sidereal  days  is  clearly  removed  by  adding  a  model  FCN  term. 
See  Figure  2  for  model  A.  From  the  analysis  of  WRMS  (Weighted  Root  Mean 
Square)  and  spectrum  and  number  of  parameters,  we  concluded  that  model  A 
(WRMS  =  0.3030  mas)  is  better  than  the  model  B  (WRMS  =  0.3051  mas).  In 
conclusion  we  adopted  the  model  of  the  FCN  term  as 

^t/’FCN  sin  Co  =  e  ^fcn* (CpcN  cos  ^FCN^  —  C'fCN  sin  ^FCnO’  (1) 

A^FCN  =  e  (x'FCNt(C'pcN  cos  WFCN^  +  Cfcn  siQ  wfcn0-  (2) 


3.  Method 

We  developed  a  new  nutation  series  by  incorporating  the  above  model  FCN 
term.  Usually  the  theory  of  the  forced  nutation  is  obtained  as  a  combination  of 
the  rigid-Earth  nutation  theory  and  the  transfer  function  of  the  nonrigid  Earth. 
As  the  former,  we  adpopted  the  RDAN98  (Roosbeek  and  Dehant,  1998).  Then 
we  removed  from  RDAN98  the  periodic  terms  whose  periods  are  longer  than  150 
years.  As  the  transfer  function,  we  selected  the  same  functional  form  as  that  of 
Herring  (1995); 


g(uj)  =  R  +  R'oj  +  „  ^FCN — 

Li>  -  ^FCN 


Rcw 

6j  -  6j cw  ’ 


(3) 
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Periodgram  of  IERS96  Periodgram  of  IERS96  +  Model  A 


Figure  2.  Spectra  of  Complex  Residuals  for  IERS96  Nutation  Theory 
with  and  without  FCN.  The  strong  peak  in  the  spectrum  of  the  IERS96 
nutation  theory  around  —400  sidereal  days  is  clearly  removed  by  adding 
a  model  FCN  term. 


where  g(uj)  is  the  ratio  of  the  nutation  amplitude  of  the  nonrigid  Earth,  <2>  is 
the  frequency  normalized  by  the  mean  angular  velocity  of  the  Earth  rotation 
Oo  as  4)pcN  =  ^fcn/^o  , where  Oq  =  7.292115  X  10_5rads_1.  Here  we  ignored 
the  contribution  of  the  inner  solid  core  of  the  Earth.  We  convolved  the  transfer 
function  with  the  nutation  theory  of  the  rigid  Earth  by  numerical  convolution 
in  the  time  domain  (Sliirai  and  Fukusliima,  2000).  We  used  the  VLBI  data 
compiled  by  USNO  for  the  period  MJD  44089.994  to  51618.250.  We  removed 
the  geodesic  nutation  discussed  by  Fukusliima  (1991)  from  the  VLBI  data  be¬ 
fore  analysis.  We  fixed  i?cw  and  u>cw  at  their  geophysically  estimated  values. 
We  also  included  in  the  parameters  to  be  fitted  the  correction  to  the  IAU1976 
precession  formula  (Lieske  :j®7  a/.,  1977),  namely  the  offsets  and  rates  of  nuta¬ 
tion.  We  included  the  complex  amplitude  of  FCN,  C'fcn  =  Cfgn  +  *^FGN  in  tlie 
parameters  to  be  fitted.  In  summary,  we  simutaneously  estimated  seven  com¬ 
plex  parameters  comprised  of  (1)  i?,  R\  i?pcNi  and  cDpcNi  (2)  two  complex  ones 
corresponding  to  the  offsets  and  rates  of  nutations  in  longitude  and  in  obliquity, 
and  (3)  the  complex  amplitude  of  the  FCN,  C'fcn-  The  estimation  was  done  by 
the  weighted  nonlinear  least-squares  method.  This  approach  is  almost  the  same 
as  that  used  to  create  the  IERS96  nutation  theory.  The  main  differences  are  in 
(1)  the  FCN  term  which  we  included  while  Herring  did  not,  (2)  the  additional 
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Periodgram  of  New  series  Periodgram  of  New  series 
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annual  progra.de  term  which  Herring  estimated  while  we  did  not,  and  (3)  the 
period  of  VLBI  observation. 


4.  Results  and  Discussion 

The  WR.MS  of  the  residuals  for  the  new  series  is  0.3075  mas,  which  is  about 
6%  smaller  than  0.325  mas  of  the  IERS96  nutation  theory.  The  estimated  pa¬ 
rameters  are  listed  in  Table  1.  Figure  1  shows  the  residuals  for  the  new  series. 
Figure  3  shows  the  periodgram  of  the  new  nutation  series  for  the  progra.de  and 
the  retrograde  bands.  The  original  R.DAN98  contains  1501  terms.  After  con¬ 
volution  with  the  transfer  function,  we  sorted  the  series  in  amplitude  and  kept 
only  the  top  180  terms.  This  is  because  the  WR.MS  of  0.3120  mas  was  only 
1.5  %  larger  than  the  0.3075  mas  value  of  the  full  series.  See  Figure  4  for  the 
number-of-terms  dependence  of  WR.MS  and  amplitude.  Clearly  the  R.MS  does 
not  siginihcantly  decrease  after  the  number  of  terms  exceeds  180.  On  the  other 
hand,  the  number-of-terms  dependence  of  the  new  series  decreases  linearly  in  a. 
log-log  plot.  In  summary,  we  consider  the  truncated  series  is  precise  enough  for 
practical  purposes.  The  full  paper  will  be  submitted  to  Astron.  -J.  soon. 
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Table  1.  Estimated  Parameters 


Real  part 

Imaginary  part 

R 

1.049018  ±  0.000011 

0.001560  ±  0.000011 

R' 

-0.2570  ±  0.0007 

-0.0406  ±  0.0007 

i?FGN 

(-1.1499  ±  0.0002)  x  10-4 

(-2.17  ±0.02)  x  10-6 

CJFCN 

(-2.3206  ±  0.0026)  X  10“3 

(3.06  ±0.26)  X  10“5 

-78  ±3 

-0.6  ±  3 

Offset  [masj 

Rate  [mas/yrj 

A '(/’ 

-42.458  ±  0.033 

-2.9297  ±  0.0047 

Ac 

-5.186  ±  0.013 

-0.2430  ±0.0019 

Number-of-Terms  ol  Dependence  of  WRMS 


Figure  4.  Number-of-terms  dependence  of  WRMS  and  nutation  am¬ 
plitude. 


Improvement  of  Nonrigid-Earth  Nutation  Theory 


229 


References 

Dehant,  V.,  et  al. ,  1999,  Celest.  Mech.  Dyn.  Astron .,  72,  245. 

Fukushima,  T.,  1991,  Astron.  Astrophys .,  244,  Lll. 

Gegout,  P.  et  al.,  1998,  Phys.  Earth.  Planet.  Inter.,  106,  337 
Herring,  T.  A.,  1995,  Hightlights  of  Astronomy,  10,  222. 

Lieske,  J.  H.,  Lederle,  T.,  Fricke,  W.,  and  Morando,  B.,  1977,  Astron.  Astro¬ 
phys.,  58,  1. 

McCarthy,  D.  D.,  1996,  HERS  Conventions’,  IERS  Technical  Note  21. 
McCarthy,  D.  D.,  2000,  private  communication. 

Roosbeek,  F.  and  Dehant,  V.,  1998,  Celest.  Mech.  Dyn.  Astron.,  70,  215. 
Sasao,  T.  and  Wahr,  J.  M.,  1981,  Geophys.  J.  R.  Astron.  Soc.,  64,  729. 
Seidelmann,  P.  K.,  1982,  Celest.  Mech.,  27,  79. 

Shirai,  T.  and  Fukushima,  T.,  2000,  Astron.  J.,  119. 

Toomre,  A.,  1974,  Geophys.  J.  R.  Astron.  Soc.,  38,  335. 


